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I. INTRODUCTION 

The idea of the possible existence of the collective 
motion in deformed nuclei similar to the scissors mo- 
tion continues to attract the attention of physicists 
who extend it to various kinds of objects, not nec- 
essary nuclei, (for example, magnetic traps, see the 
review by Heyde at al and invent new sorts of 
scissors, for example, the rotational oscillations of 
neutron skin against a proton- neutron core Q. 

The nuclear scissors mode was predicted as 
a counter-rotation of protons against neutrons in de- 
formed nuclei. However, its collectivity turned out 
to be small. From RPA results which were in quali- 
tative agreement with experiment, it was even ques- 
tioned whether this mode is collective at all @, H[ . 
Purely phenomenological models (such as, e.g., the 
two rotors model and the sum rule approach [Io| 
did not clear up the situation in this respect. Finally 
in a very recent review [l[ it is concluded that the 
scissors mode is "weakly collective, but strong on the 
single-particle scale" and further: "The weakly col- 
lective scissors mode excitation has become an ideal 
test of models - especially microscopic models - of 
nuclear vibrations. Most models are usually cali- 
brated to reproduce properties of strongly collective 
excitations (e.g. of J T = 2 + or 3~ states, giant reso- 
nances, ...). Weakly-collective phenomena, however, 
force the models to make genuine predictions and 
the fact that the transitions in question are strong 
on the single-particle scale makes it impossible to 
dismiss failures as a mere detail, especially in the 
light of the overwhelmin g ex perimental evidence for 
them in many nuclei [mjl2l]." 

The Wigner Function Moments (WFM) or phase 
space moments method turns out to be very useful 
in this situation. On the one hand it is a purely mi- 
croscopic method, because it is based on the Time 



Dependent Hartree-Fock (TDHF) equation. On the 
other hand the method works with average values 
(moments) of operators which have a direct rela- 
tion to the considered phenomenon and, thus, make 
a natural bridge with the macroscopic description. 
This makes it an ideal instrument to describe the 
basic characteristics (energies and excitation prob- 
abilities) of collective excitations such as, in par- 
ticular, the scissors mode. Our investigations have 
shown that already the minimal set of collective vari- 
ables, i.e. phase space moments up to quadratic or- 
der, is sufficient to reproduce the most important 
property of the scissors mode: its inevitable coex- 
istence with the Iso Vector Giant Quadrupole Reso- 
nance (IVGQR) implying a deformation of the Fermi 
surface. 

Further developments of the Wigner Function 
Moments method, namely, the switch from TDHF 
to Time Dependent Hartree-Fock Bogoliubov (TD- 
HFB) equations, i.e. taking into account pair corre- 
lations, allowed us to improve considerably t he q uan- 
titative description of the scissors mode [lj, [l4| : for 
rare earth nuclei the energies are reproduced with 
~ 10% accuracy and B(M1) values were reduced by 
about a factor of two with respect to their non su- 
perfluid values. However, they remain about two 
times too high with respect to experiment. We have 
suspected, that the reason of this last discrepancy is 
hidden in the spin degrees of freedom, which were so 
far ignored by the WFM method. One cannot ex- 
clude, that due to spin dependent interactions some 
part of the force of Ml transitions is shifted to the 
energy region of 5-10 MeV, where a 1 + resonance of 
spin nature is observed Q ■ 

In a recent paper [lj| the WFM method was ap- 
plied for the first time to solve the TDHF equations 
including spin dynamics. As a first step, only the 
spin-orbit interaction was included in the consider- 
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ation, as the most important one among all possible 
spin dependent interactions because it enters into 
the mean field. This allows one to understand the 
structure of necessary modifications of the method 
avoiding cumbersome calculations. The most re- 
markable result was the discovery of a new type of 
nuclear collective motion: rotational oscillations of 
"spin- up" nucleons with respect of "spin-down" nu- 
cleons (the spin scissors mode). It turns out that 
the experimentally observed group of peaks in the 
energy interval 2-4 McV corresponds very likely to 
two different types of motion: the conventional (or- 
bital) scissors mode and this new kind of mode, i.e. 
the spin scissors mode. 

Three low lying excitations of a new nature were 
found: isovector and isoscalar spin scissors and the 
excitation generated by the relative motion of the or- 
bital angular momentum and the spin of the nucleus 
(they can change their absolute values and directions 
keeping the total spin unchanged). In the frame of 
the same approach ten high lying excitations were 
also obtained: well known isoscalar and isovector Gi- 
ant Quadrupole Resonances (GQR), two resonances 
of a new nature describing isoscalar and isovector 
quadrupole vibrations of "spin- up" nucleons with re- 
spect of "spin-down" nucleons, and six resonances 
which can be interpreted as spin flip modes of vari- 
ous kinds and multipolarity. 

The obtained results are very interesting, however, 
they are only intermediate in our investigation of Ml 
modes. Our finite goal is to get reasonable agree- 
ment with experimental data for the conventional 
scissors mode, especially for its B(M1) factors which 
remain about two times too strong. Wc should keep 
in mind that only the standard spin-orbit potential 
was taken into account in the paper [15| . spin de- 
pendent residual interactions being completely ne- 
glected. 

The aim of this work is to get a qualitative under- 
standing of the influence of the spin-spin force on the 
new states analyzed in [TH, as, for instance, the spin 
scissors mode. As a matter of fact we will find that 
the spin-spin interaction does not change the gen- 
eral picture of the positions of excitations described 



in [15[. It pushes all levels up proportionally to its 
strength without changing their order. The most in- 
teresting result concerns the B(M1) values of both 
scissors modes - the spin-spin interaction strongly 
redistributes Ml strength in the favour of the spin 
scissors mode. This is a very promising fact, because 
it shows that after taking into account in addition 
pairing [l6j one may achieve agreement with exper- 
iment. 

One of the main points of the present work will, 
indeed, be that we will be able to give a tentative ex- 
planation of a recent experimental finding 17 | where 
the B(M1) values in 233 Th of the two low lying mag- 
netic states are inverted in strength in favor of the 
lowest, i.e., the spin scissors mode, when cranking 
up the spin-spin interaction. Indeed, the explana- 
tion with respect to a triaxial deformation given in 
17| yields a stronger B(M1) value for the higher ly- 
ing state, contrary to observation, as remarked by 
the authors themselves. 

The paper is organized as follows. In Sec. 2 the 
TDHF equations for the 2x2 density matrix are for- 
mulated and their Wigner transform is found. In 
Sec. 3 the model Hamiltonian is analyzed and the 
mean field generated by the spin-spin interaction is 
found. In Sec. 4 the collective variables are defined 
and the respective dynamical equations are derived. 
In Sec. 5 the results of our calculations of energies, 
B(M1) and B(E2) values are discussed. Lastly, re- 
marks and the outlook are given in the conclusion 
section. The mathematical details are concentrated 
in appendices [Aj [B] 



II. WIGNER TRANSFORMATION OF 
TDHF EQUATION WITH SPIN 

The TDHF equation in operator form reads [l6| 
ihp=[h,p\. (1) 

Let us consider its matrix form in coordinate space 
keeping all spin indices: 



ih < r, s|/5|r" 



>= 



E 



< r,s\h\r',s' >< r',s'\p\r",s" > - < r,s|p|r',.s' >< r', s'\h\r", s" > . (2) 



We do not specify the isospin indices in order to 
make the formulae more transparent. They will be 
re-introduced at the end. 

These equations will be solved by the method of 
phase space (or Wigner function) moments. To this 



end we will rewrite the expression ^ with the help 
of the Wigner transformation [l6[ . To make the for- 
mulae more readable we will not write out the coor- 
dinate dependence (r, p) of the functions. With the 



conventional notation the Wigner transform of @ can be written as 

t for s = — and i for s = 

2 * 2 



ihf^ = iHh n f n } + h n f x f n h ^ + ^ {h n f it } h it } _ # {{h n f it }} + & {{f n h it }} + 

ihf n = f n (h n ~h^) + ^{(h^ + h^)J^}-^ {{(htt-^./ti}} 

_ h n {f n f u) + ^ {h n (/ tt + f u )} + ^{{hP, (f^ f u )}} + .... 

where the functions h, f are the Wigner transforms of h, p respectively, {f,g} is the Poisson bracket of the 
functions / and g and {{/,.?}} is their double Poisson bracket; the dots stand for terms proportional to 
higher powers of h. The remaining two equations are obtained by the obvious change of arrows t^~M- 

It is useful to rewrite the above equations in terms of functions /+ = + f^, f~ = f>> — f^. By 
analogy with isoscalar / n + / p and isovector f n — f p functions one can name the functions / + and / _ as 
spin-scalar and spin-vector ones, respectively. We have: 

= ji h+ J + } + f {h~,f-} + Mh n , f"} + ih{h^, fti} + 
ihf- = l \{h\ /"} + ^{h-, /+} - 2h*f* + 2ftt4-/4.t + ^ {{h i\ fU }} _ # {{h U fit }} + 

MP = -h^f- + h-fi + jinn n + j{h + , f^} + ^{{h^ r }} ^{{h-j^}} + 

ihf" = h*f~ h-f* + |{^, /+} + |{/* + , f*} ~ ^{{h*, /-}} + ^{{h~, f*}} + (4) 

where h = h'^ ± h^. 



III. MODEL HAMILTONIAN 



The microscopic Hamiltonian of the model, harmonic oscillator with spin-orbit potential plus separable 
quadrupole-quadrupolc and spin-spin residual interactions is given by 



El + \ m u?v\ - nUSi 

2m 2 



Hqq H~ H ss 



(5) 



with 



Z N 



H=—2 y i j 

Hss = £ (-1)" | xEE S-^)S,(j) + \ 



Z N 

92- M (ri)<?2 M (rj) + Y ?2- M (ri)g 2/1 (rj) 



(6) 
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Z N 



► Jfa-r,-), (7) 



where N and Z are the numbers of neutrons and protons and are spin matrices [If 



So 



hfi o 



2 VO -1 



_S_/0 



(8) 



The quadrupole operator q 2fl = ^/l67r /5r 2 y 2 ^(^, </») can be written as the tensor product: <? 2 /x(r) = y/6{r < 
r} 2(U , where 



r. 



-i,ro,ri are cyclic coordinates [18( and C\^ Xv is a Clebsch-Gordan coefficient. 



A. Mean Field 

Let us analyze the mean field generated by this Hamiltonian. 

1. Spin-orbit Potential 
Written in cyclic coordinates, the spin-orbit part of the Hamiltonian reads 



1 



f h f h 
'°2 '-I72 



where \li 



and 



U = -hy/2Y t C^ la r v V a (9) 



h = ft(r Vi - riVo) = --7= (4 + '0 = h(r-iVi - nV-i) = 4, 

V2 



/_i = fi(r-_iV - r V_i) = — (l x - ily), 



1 

l x = —ih(yW z — zWy), l y — ~ih(zV x - xV 2 ), Z z = -iK(xW y — yV x ). (10) 



Matrix elements of in coordinate space can be obviously written as 

< r 1 ,s l \hi s \r 2 ,s 2 > = -^v(. r i){io(ri)[S Sl tS S2 t ~ S sU S S2i } 

+V2L 1 (r 1 )5 Slt S S2i - V2fi(ri)^ i; <i, 3t }5(ri - r 2 ). (11) 
The Wigner transform of (fTT|) reads [l5[ : 

/iL lS2 (r,p) =-^{io(r,p)[<W«t-^iA a i] + V^U(r,p^^ (12) 

where Z M = -i Q Ciu t i a fuPa- 



2. q-q interaction 

The contribution of H qq to the mean field potential is easily found by replacing one of the g 2 /i operators 
by the average value. We have 

V W = 6 E^ 1 )^^ ® r} 2 „. (13) 
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Hero 

= ^ + . Zfi = Kflft + Rltit) = J d(p, r){r ® r} v r+(r, p, t) (14) 

with J d(p, r) = (27r?i)~ 3 J d 3 p J d 3 r and r being the isospin index. 



3. Spin-spin interaction 



The analogous expression for H ss is found in the standard way, with the Hartree-Fock contribution given 
M by: 

Tfcfc' (*) = E Vki'k'iPw (*), (15) 

where Vki'k'i is the antisymmctrized matrix element of the two body interaction v(l,2). Identifying the 
indices k,k',l,l' with the set of coordinates (r, s,r), i.e. (position, spin, isospin), one rewrites (|15|) as 

^ F (r 1 , Sl ,r 1 ;r' 1)S ; ) r(;i) = |d 3 r 2 E E 

< ri,si,Ti;r 2 ,s 2 ,r 2 |«|r' 1 ,s / lJ Ti;r2,4)' r 2 >a.a. p( r ' 2 , 4> T 25 r 2, s 2, t 2 ; t). 
Let us consider the neutron-proton part of the spin-spin interaction. In this case 

l 

v = v(h -r 2 ) E (- 1 ) li S-^l)§ li (2)S np 5 nn , 

where fi is the position operator: f i |ri >= ri|ri >, < ri|f ilr^ >=< ri|r^ > r[ = S(ri — ri)r[. 
For the Hartree term one finds: 

< Yi,s 1 ,t 1 ;y 2 ,s 2 ,t 2 \v\y' 1 ,s' 1 ,t[]y' 2 ,s' 2 ,t! 2 >= 5(n - r[)5(r 2 - r' 2 )v(ri - r' 2 ) 
l 

E < s 1' t i; s 2,t-2|S'- a1 (1)5' ai (2)5 t1p (5 T2 „|s' 1 ,t{;s 2 ,T2 >, 



V B (r 1 ,s 1 , n ;r' 1 , S ' 1 ,r[;t) = J d 3 r 2 J d 3 r' 2 E E 

< n,si,Ti;r 2 ,s 2 ,T 2 |«|r' 1 ,s' lJ T{;r 2 ,s 2 ,r^ > p(r' 2 , s 2 , t' 2i r 2 , s 2 , r 2 ; t) 
i 

= «n P ^p E E < s iI^-m(1)I4 >< S 2 |4(2)I4 > 

<5(ri - r'J Jd 3 r 2 v(ri - r 2 ) ( o(r 2 , s' 2 , n; r 2 , s 2 , n; t). 



The Fock term reads: 



< y 1 ,s 1 ,t 1 ;y 2 ,s 2 ,t 2 \v\y' 2 ,s' 2 ,t' 2 ;y' 1 ,s' 1 ,t' 1 >= 5(y 1 - y' 2 )8{y 2 - ri)t>(4 - 
l 

E < si,Ti;s 2 ,r 2 |5_ A1 (l)4(2)^r 1 p'5r2«|s 25 ' r 2; s i 5 ' r i >, 



G 



V F (r 1 , Sl ,n;v' 1 ,s , 1 ,T , 1 ;t) = - fd 3 r 2 J d\' 2 £ £ 



< r 1 ,s 1 ,r 1 ;r2 ) S2,T2|u|r 2 ,s 2 ,T 2 ;r 1 ,s 1 ,T 1 > p(r 2 , s 2 , t' 2 ; r 2 , s 2 , r 2 ; t) 
1 

w(ri - ri)p(n,s 2 ,p; ri,s 2 ,n;t). 



Taking into account the relations 



< s|£_i|s' >= -j=6 ai 6 a >i, < s\S \s' >= ^S s ^(S st - 6 si ), < s\Si\s' >= —^=5^5*1 
and v(r — r') = x^( r ~~ r ') one finds for the mean field generated by the proton-neutron part of H ss : 



T m (r,s,T;r',s',T';t) = X~£\ 8 Tp 6 T , p 6 s i5 s >fp(r,l,n;r' ,^,n;t) + 8 s f5 s >ip(T, f, n; r',i, n;t) 



S T pS T ' n 



SsiSs>iP(r, t: P; r ', t, n; t) + <J st £ S 'tP( r > |, p; r', |, n; t) 



+-S rp 8 T < p (5 st ^ S 't - *4^s'4.) p(r, t, «; r ', t> ™; *) - /°( r J 4-1 r ', 4, *) 



S st S s/l p(r, f,p; r', |, ra; t) + <5 s4 A> t/ o(r, r', f, ra; t) 



-5 st 5 s / t p(r,t,p;r',t,n;t) - 8 s i6 s >ip(T,l,p;r',l,n;t) \6(r - r') + <p <-> n U(r - r'). 



(16) 



The expression for the mean field r pp (r, s, r; r', s', r'; i) generated by the proton-proton part of H ss can be 
obtained from Q16p by replacing index n by p and the strength constant x by x- The proton mean field is 
defined as the sum of these two terms with r = r' = p. Its Wigner transform can be written as 



V* s '(v, t) = 3 X - {5 si S s , t nf + 6 st 6 s , in f - 5 sl 6 s>l n^ - 6 st S s , t n^} 
h 2 



where n s T s (r, t) = 



d 3 p 



(17) 



-/^ s (r, p, t). The Wigner transform of the neutron mean field V£ s is obtained 



from (|17|) by the obvious change of indices pHn. The Wigner function / and density matrix p are connected 
by the relation /**'(r,p,i) = / d 3 g e~' tpq// 'V ) (ri, s, r; r 2 , s', r'; i), with q = ri - r 2 and r = |(ri + r 2 ). 
Integrating this relation over p with t' = r one finds: 

< s '(r,t) = p(r,s,T;r,s',T;t). 

By definition the diagonal elements of the density matrix describe the proper densities. Therefore n* s (r, t) 
is the density of spin- up nucleons (if s =f) or spin-down nucleons (if s =J,) . Off diagonal in spin elements of 
the density matrix n^f (r, t) are spin-flip characteristics and can be called spin-flip densities. 



IV. EQUATIONS OF MOTION 



Integrating the set of equations ((4]) over phase space with the weights 

W ={r®p}xn,{r®r}x tl ,{p®p}x l i, and 1 



(18) 



one gets dynamic equations for the following collective variables: 

LH(t) = Jd{ Vl v){r®p} Mi r{r,v,t), iZ£(t) = f d(p, r){r ® r} v /"(r, p, t), 

= fd(p,r){p®p} x „r%r,p,t), F^(t) = J d(p, r)/"(r, p, t), (19) 



where ?= +, — , t-J-i it ■ We already called the functions /+ = + and / = — spin-scalar and 
spin- vector ones, respectively. It is, therefore, natural to call the corresponding collective variables X^ (t) 
and X^~ (t) spin-scalar and spin- vector variables. The required expressions for h , tv^ and hr' are 



K = V - + mu?r 2 + 12£(-l)"SJ+(t){r ® r} 2 _„ + V+(r,t), 



h~ = -Hr,l + V-(v,t), fttJ- = ~^_ 1+ ^(r,t), /»f = + ^+(r, t), 

where according to (|17p 

7+(r,i) = -3^- X n+(r,t), F"(r,t) = 3y X n-(r,7;) + \x*n 
V^(r,t) = 3y X ^(r,t) + yxn^(r,i), V^(r,t) = 3y X nf (r,f) + yrf(r,i) 
and the neutron potentials are obtained by the obvious change of indices p <-> n. 
The integration yields: 

o 

1 



r + 



J A / / 



^A M +1 



fit 



! P+ M + 12V5^ v^7TT{^} {Z+ ® R+} Xlt 



j=o 



+ V(A-/i)(A + /i + l)L]£ +1 + V(A + M )(A-M + l)itt-! 



-?i+{r ® V} v y + + -n {r ® V} Aai V~ + r^-fr ® V} + ® V^F^ 



— P, — m w 



+ ^Ea/^+I {2A1} {^2 + ® «7 }am - iftf M^t 



y^A,i [^M-i-^ 



J=0 



<wp 4t l - 



ti 3 r [n"{r <g> V} AAi y + + n + {r ® V} Aai V~ 



-2~ / d(p, r){r ® rf v _ ^t/Ti] j 



; p^ +1 + i2^E v^JTT {«} {z+ ® pt4 }v 



J=0 



/, 2 



-^V(A-M)(A + /i + l)L^ + —77^,! 



^,oP" 



1 

3^ 



-4^,-1^ 



(20) 



~ /A[n ti HVV 1 F + |n + {r«V} v+1 ytt]^ / d(p, r){r ® p} x , l+1 [h'f^-h^r] 



= IP 



4-T 

A/i-l 



ft 2 



-^x/(A + /x)(A-/i + l)L+ M + — 77*y A ,i [<WP~ - v^^.iP^ 
-i y'd 3 r[n^{r®V} AAl _ 1 y + +n+{r®V} AM _ 1 ^] -1 f d(p, r){r ® p}^ [h^f--h~f»] 
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pit 



R 



A// 



I? 



it 



pit 



L+ -ih- 



fiR^ + v/(A- M )(A + m + 1)^+i + v / (A + /i)(A-/x+ 1)^1-1 



m 
2 



- i^ v / (A-M)(A + M + l)^ - ^ / d(p, r){r ® r}, M+1 [/T/^ - ^/-] 
^41-1 - ^|V(A + M)(A-M + l)i?t - | / d (P' r H r ® r W-i [^/- - ^1 



Z 

P+ = -2m W 2 L+ Ai + 24y5^ v ^JTT{^}{Z+®L+} AAl 



M P A " + v/(A- M )(A + M +l)P^ +1 + v/(A + M)(A-M + l)^-i 



d 3 r [{ J+ ® V} A/t F+ + {J- ® V} AA1 ^- + 2{J^ ® V} A ^ n + 2{J^ ® V} Aai ^] , 



Z 

P A " = -2m W 2 L- Ai + 24x/5^ v ^JTT{^}{Z+®L7} AAl -z^ M P 



A /A 



J=0 



d 3 r [{J" ® V} A ^+ + { J+ ® V}^"] - 2^ / d(p, r){p ® p} v [fct^t _ ^t/tl] 

2 



3=0 



d\ [{Jti ® V} V+1 F+ + {J+ ® V} V+1 ^] - ~ / d(p,v){p^p} x ^ +1 [h-f^ - h^f-j, 

2 

-ImuPL 3 *^ + 24V5^ 727 + 1(^1} {Z 2 + ® Lf - iftj V(A + A*)(A - /x + 1)P+ 



d 3 r [{jJ-t ® V} AA1 _^+ + { J+ ® V} A ^-i^ f ] 



(21) 



j^-^^p l yf c; (r,p,t) is the current. For the sake 

of simplicity the time dependence of tensors is not written out. It is easy to see that equations (|2"TT) 
are nonlinear due to quadrupole-quadrupolc and spin-spin interactions. We will solve them in the small 
amplitude approximation, by linearizing the equations. This procedure helps also to solve another problem: 
to represent the integral terms in (|2~Tj) as the linear combination of collective variables (p~9|) . that allows to 
close the whole set of equations (|2ip . The detailed analysis of the integral terms is given in the appendix IA1 

We are interested in the scissors mode with quantum number K v = 1 + . Therefore, we only need the part 
of dynamic equations with fx = 1. 



A. Linearized equations = 1), isovector, isoscalar 

Writing all variables as a sum of their equilibrium value plus a small deviation 

^ (*) = Rll + Tl Ml (t) , P Am (t) = Pll + TV (t) , L Ml (t) = L c x l + C x , (t) 



9 



and neglecting quadratic deviations, one obtains the linearized equations. Naturally one needs to know the 
equilibrium values of all variables. Evident equilibrium conditions for an axially symmetric nucleus are: 

R+ ±1 (eq)=R+ ±2 (eq)=0, B+(eq)^0. (22) 

It is obvious that all ground state properties of the system of spin up nucleons are identical to the ones of 
the system of nucleons with spin down. Therefore 

Rx^Q) = PxM) = L-^{eq) = 0. (23) 

We also will suppose 

Lt(eq) = L\\{eq) = L*(eq) = and Rf^eq) = R*(eq) = 0. (24) 

Let us recall that all variables and equilibrium quantities R~^ (eq) and Z 20 (eq) in (|21[) have isospin indices 
t = n, p. All the difference between neutron and proton systems is contained in the mean field quantities 
^20 ( e <?) an d Vt> which are different for neutrons and protons (see eq. (fT4| and (|20|) ). 



It is convenient to rewrite the dynamical equations in terms of isovector and isoscalar variables 

R\fi = R^ + Rift, P\p, = P\fi + Pxn, = L™^ 



^X^ 



R\l 



K Xu 



P: 



X/i 



pn 
'Xp, 



pP 
'Xu' 



X/.i 



L 



Xfi ' 



It also is natural to define isovector and isoscalar strength constants n\ = ^(k — R) and kq = ^(k + k) 
connected by the relation k\ = ukq Then the equations for the neutron and proton systems arc 

transformed into isovector and isoscalar ones. Supposing that all equilibrium characteristics of the proton 
system are equal to that of the neutron system one decouples isovector and isoscalar equations. This 
approximations looks rather crude, nevertheless the possible corrections to it are very small, being of the 
order ( JV ~ Z ) 2 . With the help of the above equilibrium relations one arrives at the following final set of 
equations for the isovector system: 
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where *4i, A2 are defined in appendix iBl « = — mcj 2 /(4(5oo) H3| with uj 2 = uj 2 / (l + |<5), a_i = ai = 
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The isoscalar set of equations is easily obtained from (f2"5")) by taking a — 1 and replacing x — > — x- 



V. DISCUSSION AND INTERPRETATION OF THE RESULTS 



The energies and excitation probabilities obtained by the solution of the isovector set of equations (|2"5|) are 
given in Table [I] The used spin-spin interaction is repulsive, the values of its strength constants being taken 
from the paper [2l| . where the notation x = K s /A, x — IX was introduced. The results without spin-spin 
interaction (variant I) are compared with those performed with two sets of constants K s , q (variants II, III). 
The first set of constants (variant II) was extracted by the authors of [2l[ from Skyrme forces following the 
standard procedure, the residual interaction being defined in terms of second derivatives of the Hamiltonian 
density H(p) with respect to the one-body densities p. Different variants of Skyrme forces produce different 
strength constants of spin-spin interaction. The most consistent results are obtained with SGI, SG2 [22| 
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TABLE I: Isovector energies and excitation probabilities of 1 Er. Deformation parameter 5 = 0.25, spin-orbit 
constant r\ = 0.36 MeV. Spin-spin interaction constants are: I - K s = MeV; II - K a = 92 MeV, q = —0.8; 
III - K s — 200 MeV, q = —0.5. Quantum numbers (including indices ? = +,—, fj,, 4-t) °f variables responsible 
for the generation of the present level are shown in the first column. For example: (1, 1)~ - spin scissors, (1, 1) + - 
conventional scissors, etc.. 







Eky, MeV 






fl(Ml), ii 2 N 






B(E2), B w 




i 


II 


III 


I 


II 


III 


I 


II 


III 


(1,1)" 


1.61 


2.02 


2.34 


3.54 


5.44 


7.91 


0.12 


0.36 


0.82 


(1,1) + 


2.18 


2.45 


2.76 


5.33 


4.48 


2.98 


1.02 


1.23 


1.26 


(0,0)^ 


12.80 


16.81 


20.02 


0.01 


0.01 


0.04 


0.04 


0.13 


0.52 


(2,1)" 


14.50 


18.52 


21.90 


0.01 


0.02 


0.34 


0.03 


0.13 


4.29 


(2,2)ti 


16.18 


20.61 


24.56 


0.02 


0.23 


0.03 


0.18 


3.09 


0.44 


(2,0)^ 


16.20 


22.65 


27.67 





0.03 








0.39 


0.02 


(2,1)+ 


20.59 


21.49 


22.42 


2.78 


2.19 


1.77 


35.45 


30.47 


27.43 




0.26i 


0.26i 


0.26i 


-5.4i 


-5.4i 


-5.4i 


Oi 


Oi 


Oi 



and Sk3 [23j forces. We use here the spin-spin constants extracted from Sk3 force. Another set of constants 
(variant III) was also found by the authors of [2l| phenomenologically in the calculations with a Woods- 
Saxon potential, when there is not any self-consistency between the mean field and the residual interaction. 
We tentatively will use it, because in our case also there is no self-consistency. The strength of the spin-orbit 
interaction is taken from (24[. 

One can see from Table U that the spin-spin interaction does not change the qualitative picture of the 
positions of the excitations described in [15|. It pushes all levels up proportionally to its strength (20- 
30% in the case II and 40-60% in the case III) without changing their order. The most interesting result 
concerns the relative B(M1) values of the two low lying scissors modes, namely the spin scissors (1, 1) _ 
and the conventional (orbital) scissors (1, 1) + mode. As can be noticed, the spin-spin interaction strongly 
redistributes Ml strength in the favour of the spin scissors mode. We tentatively want to link this fact to 
the recent experimental finding in isotopes of Th and Pa [13] ■ The authors have studied deuteron and 3 He- 
induced reactions on 232 Th and found in the residual nuclei 231 ' 232 > 233 Th and 232 > 233 p a "an unexpectedly 
strong integrated strength of B(M1) = 11-15 fj? N in the E 1 = 1.0 - 3.5 MeV region". The B(M1) 
force in most nuclei shows evident splitting into two Lorentzians. "Typically, the experimental splitting 
is Awjji ~ 0.7 MeV, and the ratio of the strengths between the lower and upper resonance components 
is Bl/Bjj ~ 2". (Note a misprint in that paper: it is written erroneously B2/B1 ~ 2 whereas it should 
be B1/B2 ~ 2. To avoid misunderstanding, we write here Bl instead of B\ and Bjj instead of i?2-) The 
authors have tried to explain the splitting by a 7-dcformation. To describe the observed value of Awmi 
the deformation 7 ~ 15° is required, that leads to the ratio Bl/Bjj ~ 0.7 in an obvious contradiction with 
experiment. The authors conclude that "the splitting may be due to other mechanisms". In this sense, 
we tentatively may argue as follows. On one side, theory pB| and experiment p6| give zero value of 7- 
deformation for 233 Th. On the other side, it is easy to see that our theory suggests the required mechanism. 
The calculations performed for 233 Th give Awmi ~ 0.32 MeV and B^/By ~ 1.6 for the first variant of the 
spin-spin interaction and Aljmi ~ 0.28 MeV and Bl/Bjj ~ 4.1 for second one in reasonable agreement with 
experimental values. The inclusion of pair correlations will affect our results, but one may speculate that 
the agreement between the theory and experiment will be conserved at least qualitatively. 

The energies and excitation probabilities obtained by the solution of the isoscalar set of equations (j25|) 
are displayed in the Table [TTJ The general picture of the influence of the spin-spin interaction here is quite 
close to that observed in the isovector case. The only difference is the low lying mode marked by (1, 1) + 
which is practically insensitive to the spin-spin interaction. The negligibly small negative B(M1) value of 
spin scissors appears undoubtedly due to the lack of the self consistency in our calculations. In ref [13] 
the assignment of the resonances to be of isovector type is only tentative based on the assumption that at 
such low energies there is no collective mode other than the isovector scissors mode. However, from [ijj 
one cannot exclude that also an isoscalar spin scissors mode is mixed in. From our analysis we see that the 
isoscalar spin scissors where all nucleons with spin up counter-rotate with respect the ones of spin up comes 
more or less at the same energy as the isovector scissors. So it would be very important for the future to 
pin down precisely the quantum numbers of the resonances. 
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TABLE II: The same as in Table [I] but for isoscalar excitations. 







Eis, McV 






B(M1), /4 






B(E2), B w 






I 


II 


III 


I 


II 


III 


I 


II 


III 


(1,1)" 


1.73 


2.04 


2.40 


-0.07 


-0.05 





1.12 


0.65 


0.39 


(i,i) + 


0.39 


0.37 


0.37 


0.24 


0.24 


0.24 


117.2 


117.9 


118.3 


(0,0)* 


12.83 


15.59 


18.72 











0.66 


0.31 


0.15 


(2,1)" 


14.51 


17.40 


20.65 











0.12 


0.06 


0.03 


(2,2)ti 


16.20 


19.43 


23.09 














0.07 


0.04 


(2,0)* 


16.22 


20.09 


24.80 











0.20 


0.02 


0.01 


(2,1) + 


10.28 


11.92 


13.60 











66.50 


57.78 


50.87 


(1,0)* 


0.20i 


0.20i 


0.20i 


-O.li 


-O.li 


-O.li 


30.0i 


29.8i 


30.3i 



Let us discuss in more detail the nature of the 
predicted excitations. As one sees, the generaliza- 
tion of the WFM method by including spin dynam- 
ics allowed one to reveal a variety of new types of 
nuclear collective motion involving spin degrees of 
freedom. Two isovector and two isoscalar low lying 
cigcnfrcqucncies and five isovector and five isoscalar 
high lying cigcnfrcqucncies have been found. 

Three low lying levels correspond to the excitation 
of new types of modes. For example the isovector 
level marked by (1, 1)~ describes rotational oscilla- 
tions of nuclcons with the spin projection "up" with 
respect of nucleons with the spin projection " down" , 
i.e. one can talk of a nuclear spin scissors mode. 
Having in mind that this excitation is an isovec- 
tor one, we can see that the resulting motion looks 
rather complex - proton spin scissors counter-rotates 
with respect to the neutron spin scissors. Thus the 
experimentally observed group of 1 + peaks in the in- 
terval 2-4 MeV, associated usually with the nuclear 
scissors mode, in reality consists of the excitations 
of the "spin" scissors mode together with the con- 
ventional [l[ scissors mode (the level (1, 1) + in our 
case). The isoscalar level (1,1) _ describes the real 
spin scissors mode: all spin up nucleons (protons 
together with neutrons) oscillate rotationally out of 
phase with all spin down nucleons. 

Such excitations were, undoubtedlyproduced im- 
plicitly by other methods (e.g. RPA 0, 0, EJ H3| ) , 
but they never were analysed in such terms. It is 
interesting to note, for example, that in Q the scis- 
sors mode was analyzed in so-called spin and or- 
bital components. Roughly speaking there are two 
groups of states corresponding to these two types of 
components, not completely dissimilar to our find- 
ing. Whereas the nature of the orbital, i.e. con- 
ventional scissors is quite clear, the authors did not 
analyze the character of their states which consist 
of the spin component. It can be speculated that 
those spin components just correspond to the isovec- 
tor spin scissors mode discussed in our work here. It 
would be interesting to study whether our sugges- 



tion is correct or not. This could for example be 
done in analyzing the current patterns. 

One more new low lying mode (isoscalar, marked 
by (1, 1) + ) is generated by the relative motion of the 
orbital angular momentum and spin of the nucleus 
(they can change their absolute values and directions 
keeping the total spin unchanged). 

In order to complete the picture of the low-lying 
states, it is important to discuss the state which is 
slightly imaginary. Let us first state that the nature 
of this state has nothing to do with neither spin scis- 
sors nor with conventional scissors. It can namely be 
seen from the structure of our equations that this 
state corresponds to a spin flip induced by the spin- 
orbit potential. Such a state is of purely quantal 
character and it cannot be hoped that we can accu- 
rately describe it with our WFM approach restrict- 
ing the consideration by second order moments only. 
For its correct treatment, we certainly should con- 
sider higher moments like fourth order moments, for 
instance. The spin-orbit potential is the only term 
in our theory which couples the second order mo- 
ments to the fourth order ones. As mentioned, we 
decoupled the system in neglecting the fourth or- 
der moments. Therefore, it is no surprise that this 
particular spin flip mode is not well described. Nev- 
ertheless, one may try to better understand the ori- 
gin of this mode almost at zero energy. For this, 
we make the following approximation of our diag- 
onalisation procedure to get the eight eigenvalues 
listed in Table HI We neglect in ([25|) all couplings be- 
tween the set of variables X~£ , and the set of 

variables X^,Xj^. To this end in the dynamical 
equations for we omit all terms contain- 

ing -Xvjj , X*^ and in the dynamical equations for 

X\u>X)l w e omit all terms containing X^ ,X^ . 
In such a way we get two independent sets of dy- 
namical equations. The first one (for X^ , X7) was 
already studied in [La ], where we have found that 
such approximation gives satisfactory (in compari- 
son with the exact solution) results but must be used 
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cautiously because of the problems with the angular 
momentum conservation. The second set of equa- 
tions (for XQ,X]fy splits into three independent 
subsets. Two of them were already analyzed in (l5l | 
(it turns out that these subsets can be obtained also 
in the limit 77 — > 0, which was studied there), where 
it was shown that the results of approximate calcu- 
lations are very close to that of exact calculations, 
i.e. the coupling between the respective variables 

x li> x il and x \ui x \a is ver y weak - The onl y new 



spin-vector GQR! 



subset of equations reads: 



f-it 
'-10 

pit 



.2 V 
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(26) 



of 



these equations is 



The solution 

E = i— j=t] = i 0.255 what practically coincides 

with the number of the full diagonalisation. So the 
non-zero (purely imaginary) value of this root only 
comes from the fact that z-component of orbital 
angular momentum is not conserved (only total 
spin J is conserved). However, the violation of the 
conservation of orbital angular momentum is very 
small as can be seen from the numbers. In any 
case, we see that this spin flip state has nothing 
to do with neither the spin scissors nor with the 
conventional scissors. 

Two high lying excitations of a new nature are 
found. They are marked by (2, 1) _ and following the 
paper (27j can be called spin- vector giant quadrupole 
resonances. The isovector one corresponds to the fol- 
lowing quadrupole motion: the proton system oscil- 
lates out of phase with the neutron system, whereas 
inside of each system spin up nucleons oscillate out 
of phase with spin down nucleons. The respective 
isoscalar resonance describes out of phase oscilla- 
tions of all spin up nucleons (protons together with 
neutrons) with respect of all spin down nucleons. 

Six high lying modes can be interpreted as 
spin-flip giant monopole (marked by (0,0)^) and 
quadrupole (marked by (2,0)^ and (2,2)^) reso- 
nances. 

It is a pertinent place to make following citation 
from the review by F. Osterfeld [27| : "Similar os- 
cillations to those in isospin space are also possible 
in spin space. Nucleons with spin up and spin down 
may move either in phase (spin-scalar S=0 modes) 
or out of phase (spin-vector S=l modes). The latter 
class of states is also referred to as spin excitations 
or spin-flip transitions. " On account of our results 
in this work, the latter statement that all spin ex- 
citations are of spin-flip nature should be modified. 
We predict in this paper the existence of spin ex- 
citations of non spin-flip nature - the isovector and 
isoscalar spin scissors and the isovector and isoscalar 



VI. CONCLUDING REMARKS 



The inclusion of spin-spin interaction does not 
change qualitatively the picture concerning the spec- 
trum of the spin modes found in (l5j . It pushes all 
levels up without changing their order. However, 
it strongly redistributes Ml strength between the 
conventional and spin scissors mode in the favour 
of the last one. Our calculations did not fully con- 
firm the expectations mentioned in the introduction, 
namely that essentially only the low lying part of the 
spectrum will be strongly influenced by the spin-spin 
force. Nevertheless our results turned out to be very 
useful, because they demonstrate that the common 
effect of spin-spin interaction and pair correlations 
are able to push a substantial part of the Ml force 
out of the area of the conventional scissors mode 
what is required for the reasonable agreement with 
experimental data. 

In this respect, we should mention that we did 
not include pairing in this work. Inclusion of pair- 
ing would have complicated the formalism quite a 
bit. It shall be worked out in the future. We here 
wanted to study the features of spin dynamics in a 
most transparent way staying, however, somewhat 
on the qualitative side. That is why we did not try 
to discuss in detail possible relations with experi- 
ment or to compare with the results of other the- 
ories. Nevertheless we mentioned the quite recent 
experimental work (T^j, where for the two low ly- 
ing magnetic states a stronger B(M1) transition for 
the lower state with respect to the higher one was 
found. A tentative explanation in terms of a slight 
triaxial deformation in 17j failed. However, our the- 
ory can naturally predict such a scenario with a non 
vanishing spin-spin force. It would indeed be very 
exciting, if the results of [TtJ had already discovered 
the isovector spin scissors mode. However, much 
deeper experimental and theoretical results must be 
obtained before a firm conclusion on this point is 
possible. 



In the light of the above results, the study of spin 
excitations with pairing included, will be the natural 
continuation of this work. Pairing is important for a 
quantitative description of the conventional scissors 
mode. The same is expected for the novel spin scis- 
sors mode discussed here. The effect of pairing gen- 
erally is to push up the spectrum in energy. There- 
fore, as just mentioned, it can be expected that the 
results come into better agreement with experiment. 
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Appendix A: 



d(p,r){r ®p}^[V^f^ - K iT / n ] 



V^ s being defined in (|2T)|) . It is easy to see that 
the integral I so generate moments of fourth order. 
According to the rules of the WFM method |28[ this 
integral is neglected. 



All derivations of this section will be done in the 
approximation of spherical symmetry. The inclusion 
of deformation makes the calculations more cumber- 
some without changing the final conclusions. Let us 
consider, as an example, the integral 

I h = J d(p,r){r®p}^[h^f^ - h#fK]. 

It can be divided in two parts corresponding to the 
contributions of spin-orbital and spin-spin poten- 
tials: Ih = ho + Iss, where 



Let us analyze the integral I ss (to be definite, for 
protons). In this case 

y^(r) = 3^ X n^(r) + ^rf(r), 

^(r)=3^ X nf(r) + ^x^ T (r). 

It is seen that I ss is split into four terms of identical 
structure, so it will be sufficient to analyze in detail 
only one part. For example 



I ss4 = J d(p, r){r ® p} x ^f^ = J d 3 r{r ® J n } x ^ = ]T C^ la J <Prr v j£n*, (Al) 



where J^-(r,t) = J j^^p a f^^(r,p,t). The variation of this integral reads 



SI ssi = £ C* la f d\r v [n#(eq)5j£ + J^{eq)5n^} . (A2) 

It is necessary to represent this integral in terms of the collective variables ([TO]). This problem can not 
be solved exactly, so we will use the approximation suggested in [28] and expand the density and current 
variations as a series (see appendix |B|) . 



Let us consider the second part of integral (|A2j) . With the help of formula (|B4[) we find 

V.OL 

= -£< la /dWReg)£(-l)^ (A3) 

Let us analyze at first the more simple part of this expression: 

h,i ^ -^(-l)^iV^(t) / d 3 r Y. C ti a r v JlHeq)n + = ]T(-1)%%AV (A4) 

J v,a 

We are interested in the value of fj, — 1, therefore it is necessary to analyze two possibilities: A = 1 and 
A = 2. 
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In the case A = 1, /j, = 1 we have 



(A5) 



By definition 



^p v r' (r, P) = ~ [(V„ - V>(r, s; rV 



(2nh) 



ih 



- J2 vtlMr, «)V„#(r, s') - #(r, «')V„&(r> «)], 
fe 

where k = n, I, j, m is a set of oscillator quantum numbers, t;^ are occupation numbers, and 



<lm-s(0, ' 



(A6) 



(A7) 



are single particle wave functions, X|ct( s ) — $cr,s being spin functions. Inserting (|A6[) into ()A5|) one finds 
in 1 



2 V2 



-y&OiVo-roVOyjA] 



nljm 



(A8) 



with A = m — | and A' = m + i. Remembering the definition (jlOp of the angular momentum li = 
n(?*oVi — riVo) and using the relation [l8[ f±iY/A = "F^gV ~F + l)^A±i one transforms (|A8[) into 



ifi 1 



nljm 



X " = £ <4™ / dn.+ Mr'K^-WCf^jC?* - W + A + 1) 



2/+ 1 



V2 



cZrn (r)r f , 1 R m iW-u n, 1 1 7Z ,\, n (V) 

V / nll+^m n/i+^V / nl\l— ±\m nl\l— * 



(A9) 



As it is seen, the value of this integral is determined by the difference of the wave functions of spin-orbital 
partners {vTV)^ iU+lm — ( w ^)„j|j_i| m ! which is usually very small, so we will neglect it. The only remarkable 

contribution can appear in the vicinity of the Fermi surface, where some spin-orbital partners with j ' = I + | 
and j = — can be disposed on different sides of the Fermi surface. In reality such situation happens very 
frequently, nevertheless we will not take into account this effect, because the values of the corresponding 
integrals are considerably smaller than i?2o(eg), the typical input parameter of our model. 



Let us consider now the integral 1%,! (formula (|A4j) ) for the case A = 2, /1 = 1. We have 
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r 1 4 i {eq)+r J^{eq) 



n, 



With the help of formulae (|A6I) and (|A7|) one can show by simple algebraic transformations that 

dn ri 4 l {eq) = - Jdnr j} l (eq), 
where / d£l means the integration over angles. As a result A21 = 0. 



(A10) 



(All) 
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Let us consider the second, more complicated, part of integral I2 ■ 

h.2 = -EM^^WE^ fd^jtHeq) 

,3,7 v,a •> 

= -E(- 1 ) /3+7 ^.-7W X v(/?,7)- 

0,7 



The case A = 1, /i = 1: 



1 <9n+ 



r 3r 



(A12) 



= - j E / ^ l^ n nij Wq£ H q£^_j v^i - A)(/ + a + 1) [y ;A ^ + r iA ,Y^] r/s r 7 . (Aia) 

nljm 

The angular part of this integral is 

/dfi [Ha^a + I^YJa'] r,jr 7 - E C W7 / dn M* A + F iA ,r iA ,] {r ® r} LM 



= 7= r 

V3 



1/3,17 " 



M J 



°iA,2M "r °JA',2M 



Therefore 



*ii(&7) 



/// 
6~ 



dn+{r) 3 



5 ^1/3 ^0 

2 °l^,20°i0,20 



WA 1 /-r(A' 
°/A,20 "+" ^ZA',20 



nljm v 



- l -r 5 i--0 E 1 1 - vl c lp.2o c i 



ni 



1/3 ^iO 

10,20 



IK' 



°/A,20 "+" °/A',20 



E 



21 + 1 



dr 



-iH+im^nH+l ( r ) V nl\l-^\m^nl\l-^\ ( r ) 



(A14) 



(A15) 



One sees that, exactly as in formula (|A9I) . the value of this integral is determined by the difference of the wave 
functions of spin-orbital partners (vTZ) 2 nll+ 4m — ( v '^') n i\i-i\ m near the Fermi surface, so it can be omitted 
together with In following the same arguments. 



The case A = 2, fi = 1 can be analyzed in full analogy with formulae (|A10IAllj) that allows us to take 

X'21 = 0. 



So, we have shown that the integral I2 can be approximated by zero. Let us consider now the first part of 
the integral (|A"2|) : 



J i = E / d 3 rr v n^(eq)SJ^ = E C$ la J d 3 rr v n^(eq)n+(r) E(-l) 7 A^„ 7 (t)r 7 

1/, a. v,a 7 

= E C i A M Q £(-l) 7 *£- 7 (*) / d 3 ™^( e<7 )n+(r) E ^%{r ® r} iM . (A16) 



L,M 
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This integral can be estimated in the approximation of constant density n + (r) = uq. Then 



L,M 



It is easy to show, that R^ M (eq) = 0. Let us consider, for example, the case with L = 2: 



R 2M = j rf (P: r H r ® ^}2A// it (r, P) = / dM»" ® r^M^r 1 - 
By definition 

d 3 p 



d 3 rr 2 Y 2M n^{v). (A18) 



n ss (r) 



with 0^ defined in (|A7[) . Therefore 

nljm 



IA 



^nijm J drr A Knij(r)C 3 lA ^ i 1 ^/ A , 1 _ 1 ^a/o^A/./A' — 0, 



(A17) 



(A19) 



(A20) 



where A = m — | and A' 



n, + i. The zero is obtained due to summation over to. Really, the product 
^iT 1 1 ^'w i_ i = ± ^'^ 2 ^ > 1 " (for j = i±|) does not depend on the sign of m, whereas the Clebsh-Gordan 



' 2 2 
llA 



coefficient C^j i A < changes its sign together with m. 



Summarizing, we have demonstrated that I\ + I2 — 0, hence one can neglect the contribution of the 
integrals Ih in the equations of motion. 



It is necessary to analyze also the integrals with the weight {p <8>p}v : 

i' h = [ d( P , v){p ® P } x , [h^f* h^/n] = r so + r sa 



Again we neglect the contribution of the spin-orbital part I' so , which generates fourth order moments. For 
the spin-spin contribution, we have 

I'ssi = j d(p,r){p®p} Xfl n^(r,t)fK(r,p,t) = J dhn\\{v,t)n^ {r,t), (A21) 
where Uj^(r, i) = J ^h) 3 {p ® P}*nf^( r i Pi t) i g the pressure tensor. The variation of this integral reads: 



SI' ss4 = I d 3 r 



nJ-t(eg)OTTJ(r, t) + U^(eq)Sn^(r, t) 



t+ ( Br ,\x„Vti 



(A22) 



The pressure tensor variation is defined in appendix |Bj With formula (|B6p one finds for the first part of 
(|A22jl : 



I[ = J d 3 rn^(eq)5U\\ i (r,t) ~ T^(i) J d 3 rn it (eq)n + (r) ~ T^(t)n a j d 3 rn^ (eq) = 0. (A23) 
The last equality follows obviously from the definition of n*' (|A19|) . 



18 

The second part of (|A22[) reads: 

I' 2 = fd 3 rIlK(eq)6n^(r,t) 

= -Et- 1 )^ 3 ^ 6 ^ (A24) 

Let us consider at first the simpler part of this integral 

- Y,(-lfNg_p(t) J d\li\\{e q )n + {r). (A25) 





■ 



The value of the last integral is determined by the angular structure of the function nj^(r). We are interested 
in A = 2, fi = 1. By definition 

l4i(r)=/ j^{P®PhifHr,p)=^Ctl. la J j^ P ^fHr,p) 
= / (^3P^o/ n (r, P ) = [(Vi - Vi)(V - V„)p(r' t,r i)] r , =r 

= -^E^I V > V «^( r 't)]^(4) - [Vi^(r,t)][V ^(r,|)] 

-[Vo^^t^tVi^^i^+^^^lViVo^Cr^)]} (A26) 
with being defined by (|A7|) ■ Taking into account formulae [HI 



(Z±A+1)(Z±A + 2) Z (I T A-1)(ItA) I + 1 

V±lFiA - 2(2/ + l)(2Z + 3) ^ +1 < A±1 "V 2(2Z-1)(2/ + 1) r i " 1 ' A±1 ' 



one finds that 



/ (Z + l) 2 -A 2 Z Z 2 - A 2 Z + l 

VoFzA - "V (2TT1K2T+3) r +1 ' A + V (2Z-1)(2Z + 1) r 



/ ( eg ) n +(r) = Z* 2 £ v\ ijm J *n+(r)^,(r)(^: - tf J.*-*) (21 + 3)(2i + 1)(2I -1) m = ° (A2?) 

due to summation over m. The more complicated part of the integral (|A24[) is calculated in a similar way 
with the same result, hence I' 2 = 0. 

So, we have shown that I[ + 1' 2 — 0, therefore one can neglect by the contribution of integrals I' h (together 
with Ih) into equations of motion. 

• And finally, just a few words about the integrals with the weight {r (g) r} A/i : 

III = J d(p, r){r (» r} AM [h^f* h*fU] = C + I' s ' s . 
The spin-orbital part J" Q is neglected and for the spin-spin part we have 

I'La = J d(p, r){r <g> r}x^{r, t)/ n (r, p, t) = J d 3 r{r <g> r} A/ ^(r, t)n^(r, t). (A28) 
The variation of this integral reads: 

SI "s4 = [ d 3 r{r ® r} AA1 [n^(eg)(5n n (r,t) + n n (eq)(5n^(r, i)]. (A29) 
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With the help of formulae (|A19[) and (|B4[) the subsequent analysis becomes quite similar to that of the 
integral (|A16[) with the same result, i.e. 1% ~ 0. 

• The integrals jd[p,r)W\ ll [h~f^ — h^f~] and Jd(p,r)Wx f _ t [h~f^ — h^f~] , where is any of 
the above mentioned weights, can be analyzed in an analogous way with the same result. 



Appendix B: 

According to the approximation suggested in |28j . the variations of density, current, and pressure tensor 
are expanded as the series 



5tf(r,t) = -^(-lfV„l+(r) 

I 



^w + Et- 1 ) 7 ^^ 



(Bl) 



6JUr,t)=n + (r) 



(B2) 



<5II^(r,t)=n + (r) 



7 A', ft' 



(B3) 



Putting these series into the integrals (|A2| IA22|) . one discovers immediately that all terms containing ex- 
pansion coefficients N, K, T with odd numbers of indices disappear due to axial symmetry. Furthermore, 
we truncate these series omitting all terms generating higher (than second) order moments. So, finally the 
following expressions are used: 



«Jn s (r,t) ~ -^(-l)^V_Jn+(r)^(-l)^ !7 (t)r_ 7 ) 

= -em)" +T,(- i rNi,(t) 1 - 

P l 7 



9r 



-P r --y 



(B4) 



*J|(r,t)~n + (r)X;(-l) 7 ^,- 7 (*> 



(B5) 



and 



SUl^v,t)^n+(v)T^(t). 



(B6) 



The coefficients iVg (t) and -fT| are connected by the linear relations with collective variables 'R-\ u (t) 



V P -7 



A// ^ 



and respectively 
where 



^ = V2 fl§§ - i$g = ^§ (l + ^) - -A2 = fl^/V2 + IQ = -9^ (l - 2 -6 ) 
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Rio = Q20/V&, Roo = -Qoo/VZ, Q20 = UQq , Qqo = A< r 2 >= §ARg. 



N l-i 



AT S — 

ly -l,0 - A A > JV -1,1 



AA 2 



4Ai 



2Ai 

44i ' 



1.1 



. ^0,1 

2^2 



2A 2 
4^l 2 ' 



(B9) 



A, /i 



d 3 r{r (g> 5J ? }a m 



A^io^,„ - -42 + C^_ 1A1 Kl_^)] . (BIO) 



V3C 



2-2 ^ Vsqcu + ^-i) ^ s v^Ao + £ 20 + V2C 

A -1.0 — ' J -1,1 — 



A 2 
V2-4 2 



V^-^2,0 — -^0,0 



. ^0,1 



A! 

V3(Ai-Ai 



V^-^10 -^20 V^-^-OQ 7^; 

) "1,0 



V2^2 

^(Ai + Ai) 
V3A 



\/2^2 



\/2.4.i 



, «!,i = - 



-A 2 



(Bll) 



The coefficient T^(t) is connected with V\At) by the relation V^^t) = AT^ (t), A being the number of 
nucleons. I 
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